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SUMMARY  This  paper  proposes  a  high-speed  architecture  to  realize 
two- variable  numeric  functions.  It  represents  the  given  function  as  an  edge¬ 
valued  multiple- valued  decision  diagram  (EVMDD),  and  shows  a  system¬ 
atic  design  method  based  on  the  EVMDD.  To  achieve  a  design,  we  char¬ 
acterize  a  numeric  function  /  by  the  values  of  /  and  p  for  which  /  is  an 
/-restricted  Mp-monotone  increasing  function.  Here,  /  is  a  measure  of  sub¬ 
functions  of  /  and  p  is  a  measure  of  the  rate  at  which  /  increases  with  an 
increase  in  the  dependent  variable.  For  the  special  case  of  an  EVMDD, 
the  EVBDD,  we  show  an  upper  bound  on  the  number  of  nodes  needed 
to  realize  an  /-restricted  Mp-monotone  increasing  function.  Experimen¬ 
tal  results  show  that  all  of  the  two-variable  numeric  functions  considered 
in  this  paper  can  be  converted  into  an  /-restricted  Mp-monotone  increas¬ 
ing  function  with  p  =  1  or  3.  Thus,  they  can  be  compactly  realized  by 
EVBDDs.  Since  EVMDDs  have  shorter  paths  and  smaller  memory  size 
than  EVBDDs,  EVMDDs  can  produce  fast  and  compact  NFGs. 
key  words:  two-variable  numeric  function  generators  (NFGs),  edge- 
valued  multiple-valued  decision  diagrams  (EVMDDs),  edge-valued  bi¬ 
nary  decision  diagrams  (EVBDDs),  graph-based  representation  of  numeric 
functions,  programmable  memory-based  architecture 

1.  Introduction 

Numeric  functions  have  wide  applications  including  com¬ 
puter  graphics,  direct  digital  frequency  synthesizers  [5], 
and  digital  signal  processing.  Various  design  methods 
for  numeric  function  generators  (NFGs)  have  been  de¬ 
veloped  [18].  However,  most  existing  methods  are  in¬ 
tended  for  one-variable  numeric  functions  [7],  [16],  [21], 
[25],  [29]— [3 1  ],  and  only  a  few  methods  have  been  reported 
for  specific  multi-variable  numeric  functions  [9],  [10],  [34], 
Thus,  different  numeric  functions  require  different  methods. 
As  far  as  we  know,  no  systematic  design  method  for  generic 
multi-variable  numeric  functions  has  been  presented. 

A  straightforward  design  method  for  an  arbitrary  multi- 
variable  function  is  to  use  a  single  memory.  This  method 
produces  a  fast  NFG,  but  requires  a  2/c"-word  memory  to 
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realize  a  ^-variable  function  with  n  hits  for  each  variable. 
Thus,  even  for  a  computation  with  a  small  number  of  bits, 
this  method  is  impractical  because  of  large  memory  size. 

To  produce  practical  NFGs,  we  consider  a  design 
method  using  decision  diagrams  (DDs)  for  two-variable 
NFGs.  DDs,  such  as  binary  DDs  (BDDs),  can  compactly 
realize  various  functions  [8],  [  17],  [35]-[37].  However,  DDs 
are  not  able  to  represent  all  classes  of  the  functions  com¬ 
pactly.  Thus,  choosing  a  DD  appropriate  to  a  given  class 
of  functions  is  important.  Although  DDs  suitable  for  one- 
variable  numeric  functions  have  been  presented  [21],  [23], 
[29],  [33],  as  far  as  we  know,  no  study  on  graph-based  rep¬ 
resentations  for  multi-variable  numeric  functions  has  been 
reported. 

First,  we  present  a  DD  appropriate  for  two-variable  nu¬ 
meric  functions.  And  then,  we  propose  a  design  method 
and  an  architecture  for  two- variable  NFGs  using  this  DD. 
To  analyze  complexities  for  two-variable  numeric  functions, 
we  introduce  a  new  class  of  integer-valued  functions,  l- 
restricted  M/z-monotone  increasing  functions.  We  derive 
an  upper  bound  on  the  number  of  nodes  in  an  edge-valued 
BDD  (EVBDD)  for  this  type  of  function.  Theoretical  analy¬ 
sis  and  experimental  results  show  that  edge-valued  multiple¬ 
valued  DDs  (EVMDDs)  can  compactly  represent  both  one- 
and  two-variable  numeric  functions,  and  our  NFGs  using 
EVMDDs  can  compactly  realize  such  functions  with  the 
same  architecture. 

This  paper  is  organized  as  follows:  Section  2  intro¬ 
duces  a  fixed-point  representation  to  convert  a  real-valued 
numeric  function  into  an  integer-valued  function.  Section  3 
considers  representations  of  two-variable  numeric  functions 
using  DDs.  It  introduces  the  /-restricted  M /umonotone  in¬ 
creasing  function,  and  derives  an  upper  bound  on  the  num¬ 
ber  of  nodes  in  an  EVBDD  for  this  type  of  function.  Sec¬ 
tion  4  presents  an  architecture  and  a  design  method  for 
NFGs  based  on  EVMDDs.  Experimental  results  using  an 
FPGA  are  also  presented.  Section  5  concludes  the  paper. 

2.  Preliminaries 

2. 1  Number  Representation  and  Precision 

Definition  1:  Let  B  -  [0, 1),  Z  be  the  set  of  the  integers, 
and  R  be  the  set  of  the  real  numbers.  A  logic  function  is 
B"  — >  Bm .  An  integer- valued  function  is  B"  —>  Z.  A  one- 
variable  numeric  (real)  function  is  K.  — »  R.  And,  a  two- 
variable  numeric  function  is  R  x  R  — >  R. 
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Definition  2:  A  numeric  function  generator  (NFG)  is  a 
logic  circuit  that  computes  approximate  values  for  a  numeric 
function.  A  one-variable  NFG  is  a  logic  circuit  for  a  one- 
variable  numeric  function  fix),  whose  input  is  x,  and  out¬ 
put  is  an  approximate  value  for  fix).  A  two- variable  NFG 
is  a  logic  circuit  for  a  two-variable  numeric  function  fix,  y ), 
whose  inputs  are  x  and  y,  and  output  is  an  approximate  value 
for  fix,  y). 

Definition  3:  The  binary  fixed-point  representation  of  a 
number  is  denoted  by 

A  —  (xn_int-  |  Xn_int-2  ■  •  •  X]  A'(j ■  X  _  I  A'  -2  .  .  .  X-n_frac) 2, 

where  .a,  e  {0, 1),  nJnt  is  the  number  of  bits  for  the  inte¬ 
ger  part,  and  n_frac  is  the  number  of  bits  for  the  fractional 
part  of  X.  This  is  the  two’s  complement  representation: 

y  _  _r)ri-int—l  v  .  2  v  ,  ,  r)—ri-frac  v 

^  1  '  ^  -A-ti-int— 2  '  •  •  •  '  ^  ^-n^frac' 

To  distinguish  the  set  of  binary  variables  from  the  numeric 
value  X,  we  use  brackets  {  }.  Specifically,  {X}  denotes  the 
set  of  binary  variables  x 

Definition  4:  Precision  is  the  total  number  of  bits  in  a  bi¬ 
nary  fixed-point  representation.  Specially,  n-bit  precision 
specifies  that  n  bits  are  used  to  represent  the  number;  that  is, 
n  =  nJnt  +  n  frac.  In  this  paper,  an  n-bit  precision  func¬ 
tion  fiX,  Y)  means  that  both  of  the  input  variables  X  and 
Y  have  n-bit  precision.  The  number  of  fractional  bits  for 
function  values  is  n.  This  is  because  the  range  of  a  function 
can  be  different  than  its  domain.  Thus,  function  values  have 
inJnt  +  n)- bit  precision,  where  nJnt  is  the  number  of  integer 
bits  for  function  values. 

We  can  convert  an  n-bit  precision  two-variable  numeric 
function  into  a  2n-input  m-output  logic  function,  where 
m  =  nJnt  +  n.  The  logic  function  can  be  converted  into 
an  integer-valued  function  by  considering  binary  vectors  as 
integers.  That  is,  we  can  convert  an  n-bit  precision  two- 
variable  numeric  function  into  an  integer-valued  function: 

if"  — >  P,„,  where  Pm  =  {0, 1 _ ,2"'  -  1}.  In  this  paper, 

two-variable  numeric  functions  are  converted  into  integer¬ 
valued  functions  in  this  way.  And,  in  the  text  that  follows, 
we  assume  that  aq  and  ytt  denote  the  least  significant  bits  in 
the  fixed-point  representations  of  X  and  Y,  respectively. 

Example  1:  Table  1  (a)  is  the  function  table  for  the  Eu¬ 
clidean  norm  function  fX2  +  Y2  for  a  two-dimensional  vec¬ 
tor  from  (0, 0)  to  (X,  Y).  The  2-bit  precision  fixed-point  rep¬ 
resentation  of  this  function  is  the  logic  function  fifX,  Y)  in 
Table  1  (b).  By  converting  output  vectors  into  integers,  we 
have  the  integer- valued  function  fiX,  Y)  of  fbiX,  Y)  in  Ta¬ 
ble  1  (c).  In  this  paper,  the  2-bit  precision  2-D  norm  function 
denotes  the  integer-valued  function  f{X,  Y )  in  Table  1  (c). 

(End  of  Example) 

2.2  Decision  Diagrams 

This  subsection  summarizes  the  DDs  used  in  this  paper.  For 
more  detail  on  definitions  and  reduction  rules,  see  [8],  [27], 
[37], 


Table  1  Tables  for  2-bit  precision  2-D  norm  function. 

(a)  Table  for  2-D  norm,  (b)  Table  for  fb(X,  Y ).  (c)  Table  for  f(X,  Y ). 
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Definition  5:  A  multi-terminal  binary  decision  diagram 
(MTBDD)  [6]  is  an  extension  of  a  BDD  [3],  [17],  and  rep¬ 
resents  an  integer-valued  function.  In  the  MTBDD,  the  ter¬ 
minal  nodes  are  labeled  by  integers. 

Definition  6:  A  binary  moment  diagram  (BMD)  [4]  is 
a  rooted  directed  acyclic  graph  (DAG)  representing  an 
integer- valued  function.  The  BMD  is  obtained  by  repeatedly 
applying  the  arithmetic  transform  expansion  f  -  fo  +  x,(/i  ~ 
/o)  to  the  integer-valued  function,  where  fo  =  /(x,-  =  0),  and 
fi  =  fix i  =1).  The  BMD  consists  of  terminal  nodes  repre¬ 
senting  the  arithmetic  coefficients,  and  non-terminal  nodes 
representing  the  arithmetic  transform  expansions.  Each  non¬ 
terminal  node  has  two  edges  corresponding  to  two  terms:  fo 
and  Xiifi  -  fo )  in  the  arithmetic  transform  expansion. 

Definition  7:  An  edge-valued  BDD  (EVBDD)  [14],  [15]  is 
a  variant  of  a  BDD,  and  represents  an  integer-valued  func¬ 
tion.  The  EVBDD  is  obtained  by  repeatedly  applying  the 
expansion  /  =  x//o  +  xff[  +  a)  to  the  integer-valued  func¬ 
tion,  where  fo  =  fixt  =  0),  f[+a-f\  =  /(.A;  =  1),  and  a 
is  the  constant  term  of  f\ .  The  EVBDD  consists  of  only  one 
terminal  node  representing  0  and  non-terminal  nodes  with 
1 -edges  having  integer  weights  a.  In  the  EVBDD,  0-edges 
always  have  zero  weights.  The  incoming  edge  to  the  root 
node  can  have  a  non-zero  weight.  In  a  reduced  EVBDD, 
each  node  represents  a  distinct  sub-function. 

Definition  8:  For  an  n-bit  precision  number  X,  if  ]A]  = 
{X„}U{XH_1}U. .  .U{Zi],  {Xf  *  0,  and{X;]npr7]  =  0  (i  *  j), 
then  [Xu,  Xu-\ , . . . ,  X\]  is  a  partition  of  X.  Each  X,  forms  a 
super  variable.  Let  |A,|  =  k,  and  ku  +  ku^\  +  . . .  +  k\  =  n. 
Then,  by  considering  each  super  variable  as  a  multiple¬ 
valued  variable,  an  integer- valued  function  fiX)  :  B"  — >  Z 
can  be  converted  into  a  multiple-valued  input  integer  func¬ 
tion  fiXu,  A„_| , . . . ,  X\)  :  Pu  X  Pu  \  x  ...  x  Pi  — >  Z,  where 
Pi  =  {0, 1,2,..., 2*'  -  1). 

Definition  9:  An  edge-valued  multiple-valued  decision  di¬ 
agram  (EVMDD)  [21]  is  an  extension  of  an  MDD  [13],  and 
represents  a  multiple-valued  input  integer  function.  It  con¬ 
sists  of  one  terminal  node  representing  0  and  non-terminal 
nodes  with  edges  having  integer  weights,  and  0-edges  al¬ 
ways  have  zero  weights.  As  shown  in  Fig.  1,  an  EVMDD  is 
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Fig- 1  Conversion  of  EVBDD  nodes  into  an  EVMDD  node. 


(c)  EVBDD 


Fig.  2  Four  types  of  DDs  for  the  2-bit  precision  2-D  norm  function. 


obtained  by  merging  non-terminal  nodes  in  an  EVBDD  ac¬ 
cording  to  the  partition  of  X.  In  an  EVMDD,  each  multiple¬ 
valued  variable  may  have  a  different  domain. 


3.  Graph-Based  Representations  of  Two- Variable  Nu¬ 
meric  Functions 

This  section  introduces  an  /-restricted  Mp-monotone  in¬ 
creasing  function,  and  derives  an  upper  bound  on  the  num¬ 
ber  of  nodes  in  an  EVBDD  for  the  /-restricted  Mp-monotone 
increasing  function.  Experimental  results  in  this  section 
show  that  EVBDDs  for  two-variable  numeric  functions  are 
more  compact  than  MTBDDs  and  BMDs. 

3.1  /-Restricted  Mp-Monotone  Increasing  Functions 

Definition  10:  An  n-bit  precision  integer-valued  function 
/(X)  such  that  0  <  f(X  +  1)  -  f(X)  <  p  and  /( 0)  =  0  is 
a  totally  Mp-monotone  increasing  function  (or  simply,  Mp- 
monotone  increasing  function).  Here,  X  +  1  is  the  binary 
representation  of  the  independent  variable.  That  is,  for  an 
Mp-monotone  increasing  function  f(X),  f( 0)  =  0,  and  the 
increment  of  X  by  one  increases  the  value  of  /(X)  by  at  most 
72- 

Adding  1  as  in  X  +  1  is  simply  incrementing  the  stan¬ 
dard  binary  number  X.  It  should  not  be  confused  with 
adding  1  to  a  real- valued  variable  x. 

Definition  11:  An  «-bit  precision  integer-valued  function 
/(X)  is  an  /-restricted  Mp-monotone  increasing  function 
when,  for  1  <  /  <  n,  all  the  /-bit  precision  sub-functions 
g(Xj)  of  /  are  Mp-monotone  increasing  functions,  where 
{X)  =  {xn-ux„-2,...,x0},  (X,)  =  {x/_i,x/_2,  ...,x0],  and 
g(X/)  =  f{a,  X/),  for  all  assignments  a  of  values  to 

(Xti—  1  Xji — 2  -  -  -  A/  )o . 

Theorem  1:  For  an  77-bit  precision  /-restricted  Mp- 
monotone  increasing  function  /(X),  the  number  of  nodes  in 
the  EVBDD  is  at  most 


Example  2:  Figure  2  (a),  (b),  (c),  and  (d)  show  the 
MTBDD,  the  BMD,  the  EVBDD,  and  the  EVMDD  of  the 
2-bit  precision  2-D  norm  function  in  Table  1  (c).  For  read¬ 
ability  of  the  figures,  several  terminal  nodes  are  not  shared. 
In  Fig.  2  (a)  and  (c),  dashed  lines  and  solid  lines  denote 
0-edges  and  1-edges,  respectively.  Note  that  the  EVBDD 
has  weighted  1-edges.  In  Fig.  2(b),  ‘A’  in  a  circle  denotes 
the  arithmetic  transform  expansion.  And,  in  Fig.  2  (d),  the 
set  of  binary  variables  {X)  U  {Y}  is  partitioned  into  {X2)  = 
{xi,x0,yi}  and  {Xi }  =  {yol-  To  obtain  the  function  value 
3  for  X  =  (10)2  and  Y  =  (10)2,  in  the  MTBDD,  we  tra¬ 
verse  the  MTBDD  from  the  root  node  to  a  terminal  node 
according  to  the  input  values,  and  obtain  the  function  value 
from  the  terminal  node.  In  the  BMD,  we  obtain  the  func¬ 
tion  value  by  computing  the  arithmetic  transform  expan¬ 
sion  f  =  fo  +  Xj(f\  -  fo)  recursively  at  each  non-terminal 
node.  And,  in  the  EVBDD  and  the  EVMDD,  we  obtain 
the  function  value  as  the  sum  of  the  weights  for  the  edges 
traversed  from  the  root  node  to  the  terminal  node.  Note 
that  we  traverse  the  EVMDD  using  X2  =  5  and  Xi  =  0. 

(End  of  Example) 


1 

2n-‘  +  J](p+lf-1  -L  (1) 

7=1 

where  /  is  the  largest  integer  satisfying  2"~l  >  (p+l)2*-1,  and 
the  variable  order  of  the  EVBDD  is  x„-i,x„-2, . . . ,  xo  (from 
the  root  node  to  the  terminal  node). 

Proof:  See  the  appendix. 

Theorem  1  also  holds  for  EVMDDs  because  an 
EVMDD  is  obtained  by  merging  non-terminal  nodes  in 
an  EVBDD.  Note  that  the  upper  bound  for  /-restricted 
Mp-monotone  increasing  functions  shown  in  Theorem  1 
is  equal  to  the  upper  bound  for  totally  Mp-monotone  in¬ 
creasing  functions  shown  in  [21].  This  upper  bound  is 
much  smaller  than  the  worst-case  upper  bound,  2",  which 
is  reached  by  EVBDDs  for  power  functions  and  polynomial 
functions  [23]. 

Example  3:  Consider  a  16-bit  precision  /-restricted  Mp- 
monotone  increasing  function.  If  p  =  1,  then  we  have  /  =  3, 
and  the  upper  bound  given  by  (1)  is  8,327.  If  p  =  3,  then 
1  =  2,  and  the  upper  bound  is  16,450.  (End  of  Example) 
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Definition  12:  An  n-bit  precision  integer-valued  function 
/(X)  is  an  extended  /-restricted  Mp-monotone  increasing 
function  when,  for  1  <  /  <  n,  all  the  /-bit  precision 
sub-functions  of  /  are  Mp-monotone  increasing  functions 
g(X/ )  or  represented  by  g{X{)  +  b ,  where  b  is  an  inte¬ 
ger,  {X}  =  {xn-i,Xn-2,...,X0},  {X,}  =  {X/_1,A/_2,...,X0}, 
and  sub-functions  are  f(ct,  X/),  for  all  assignments  a  to 
(xn—\  Xn—2  .  .  .  A7  )2  • 

Lemma  1:  Let  /(X)  be  an  extended  /-restricted  Mp- 
monotone  increasing  function.  For  any  integer  /'  satisfying 
1  </'</,  /(X)  is  an  extended  /'-restricted  Mp-monotone 
increasing  function. 

Proof:  This  follows  from  Definition  12.  I 

Lemma  2:  Let  /(X)  be  an  /-restricted  Mp-monotone  in¬ 
creasing  function,  and  let  g(X)  be  an  extended  /-restricted 
Mp-monotone  increasing  function  that  is  obtained  by 
adding  constant  values  to  the  /-bit  precision  sub-functions 
of  f.  Then,  the  EVBDDs  for  /(X)  and  g(X)  have  the  same 
number  of  nodes. 

Proof:  See  the  appendix. 

Corollary  1:  Let  /(X)  be  an  extended  /-restricted  M p- 
monotone  increasing  function,  and  let  g(X)  be  an  affine 
transformation  off:  g(X)  =  af(X)  +  b,  where  a  and  b  are  in¬ 
tegers.  Then,  the  EVBDDs  for  f(X)  and  g(X)  have  the  same 
number  of  nodes. 

EVBDDs  can  compactly  represent  not  only  /-restricted 
Mp-monotone  increasing  functions,  but  also  their  extended 
classes  of  functions.  This  property  is  helpful  to  compactly 
represent  various  two- variable  numeric  functions. 

3.2  Two- Variable  Numeric  Functions 

As  shown  in  Sect.  2,  n-bit  precision  two-variable  numeric 
functions  can  be  converted  into  2n-bit  precision  integer¬ 
valued  functions.  That  is,  n-bit  precision  two-variable  func¬ 
tions  /(X,  T)  can  be  converted  into  2n-bit  precision  one- 
variable  functions  /(Z),  where 

Z  =  2"X  +Y-  (x„_i  xn-2  ...  xq  v„-i  y„-2  ■  ■ .  v0)2. 

When  /(Z)  is  an  /-restricted  Mp-monotone  increasing  func¬ 
tion  for  the  largest  integer  /  satisfying  22"_/  >  (p  +  1)2~\ 
Theorem  1  gives  the  upper  bound  on  the  number  of  nodes  in 
an  EVBDD  for  /(X,  Y). 

Example  4:  As  shown  in  Table  2  (a),  the  2-bit  precision 
2-D  norm  function  Vx2  +  Y2  can  be  converted  into  an  ex¬ 
tended  2-restricted  Ml-monotone  increasing  function  f(Z). 

Table  2  Function  tables  for  2-bit  precision  two-variable  functions. 

(a)  Table  for  2-D  norm,  (b)  Table  for  yyy .  (c)  Table  for  g(Z). 
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3 

4 

2 

0 

t 

1 

2 

to 

0 

-1 

-1 

-2 

3 

3 

3 

4 

4 

3 

0 

t 

1 

2 

11 

0 

-1 

-1 

-2 

Note  that,  in  the  table,  values  increase  by  at  most  one  for 
each  column.  Similarly,  the  2-bit  precision  two-variable 
function  yy  shown  in  Table  2(b)  can  be  converted  into 
an  affine  transformation  of  the  extended  2-restricted  Ml- 
monotone  increasing  function  g(Z)  shown  in  Table  2(c): 
-lx  g(Z).  (End  of  Example) 

We  now  show  two-variable  numeric  functions  that  can 
be  converted  into  integer-valued  functions  of  classes  dis¬ 
cussed  in  Sect.  3.1.  Asa  result,  their  EVBDDs  are  small. 

Lemma  3:  Let  h(  Y)  be  an  n-bit  precision  Mp-monotone 
increasing  function.  Then,  for  an  arbitrary  one-variable 
function  g(X),  a  two-variable  function  /(X,  T)  =  g(X)  +  h(  Y) 
is  an  extended  n-restricted  Mp-monotone  increasing  func¬ 
tion. 

Proof:  See  the  appendix. 

Lemma  4:  Let  h{Y)  be  an  n-bit  precision  Mp-monotone 
increasing  function.  Then,  for  an  arbitrary  one-variable 
function  g(X),  the  two-variable  function  /(X,  Y)  =  g(X)  - 
h(Y )  can  be  converted  into  an  affine  transformation  of  an 
extended  n-restricted  Mp-monotone  increasing  function. 
Proof:  See  the  appendix. 

Lemma  5:  Let  h(  Y)  be  an  n-bit  precision  Mp-monotone 
increasing  function,  and  let  g(X)  be  a  real  function  satisfying 
0  <  g(X)  <  1.  Then,  an  n-bit  precision  two-variable  func¬ 
tion  /(X,  Y)  =  g(X )  •  h(Y)  is  an  n-restricted  Mp-monotone 
increasing  function. 

Proof:  See  the  appendix. 

In  Lemma  5,  if  the  range  of  g(X)  is  large,  then  the 
EVBDD  can  be  large.  For  example,  the  n-bit  multiplier  re¬ 
quires  0(2")  nodes  [36]. 

Lemma  6:  Let  h(Y)  be  an  affine  transformation  of  an  n-bit 
precision  Mp-monotone  increasing  function,  and  let  g(X) 
be  a  real  function  satisfying  0  <  g(X)  <  1.  Then,  an  n-bit 
precision  two-variable  function  /(X,  Y)  =  g(X)  ■  h(Y)  can 
be  converted  into  an  affine  transformation  of  an  extended 
n-restricted  Mp-monotone  increasing  function. 

Proof:  See  the  appendix. 

Example  5:  2-bit  precision  function  Wy  is  an  affine  trans¬ 
formation  of  an  Ml-monotone  increasing  function  [21].  As 
shown  in  Example  4,  /(X,  T)  =  yy  can  be  converted  into 
an  affine  transformation  of  an  extended  2-restricted  Ml- 
monotone  increasing  function.  (End  of  Example) 

Since  many  common  one-variable  numeric  functions 
can  be  converted  into  Mp-monotone  increasing  func¬ 
tions  [23],  many  two-variable  numeric  functions  obtained 
by  four  arithmetic  operations  of  them  can  be  converted  into 
extended  /-restricted  Mp-monotone  increasing  functions  as 
shown  in  the  above  lemmas. 

In  the  following,  we  show  that  various  two-variable 
numeric  functions,  as  well  as  the  above  functions,  can  be 
converted  into  extended  /-restricted  Mp-monotone  increas¬ 
ing  functions.  Table  3  compares  the  numbers  of  nodes  in 
MTBDDs,  BMDs,  and  EVBDDs  for  certain  8-bit  precision 
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Table  3  Numbers  of  nodes  in  MTBDDs,  BMDs,  and  EVBDDs  for  8-bit 
precision  two-variable  numeric  functions. 


Numeric 

functions 

Type  of 
function 

Number  of  nodes 

Ri 

r2 

MTBDD 

BMD 

EVBDD 

six2  +  Y2 

Ml 

12,969 

25,084 

2,566 

20 

10 

arctan(  yyy ) 

Ml+ 

8,997 

26,158 

3,134 

35 

12 

ln(X+  l)sin(y) 

Ml 

9,776 

25,994 

3,444 

35 

13 

VI  sin  (Y) 

Ml 

11,543 

26,542 

3,483 

30 

13 

sin(  six 2  +  Y2) 

Ml 

11,521 

27,858 

4,013 

35 

14 

sin  (XY) 

Ml 

11,282 

21,746 

3,789 

34 

17 

X/(Y+  1) 

Ml+ 

9,664 

25,878 

3,162 

33 

12 

XY i  six 2  +  Y2 

Ml 

9,325 

23,634 

2,269 

24 

10 

WaveRings 

M3+ 

17,423 

27,691 

5,047 

29 

18 

Average 

11,389 

25,621 

3,434 

30 

13 

Domain  of  the  functions  is  0  <  X  <  1  and  0  <  Y  <  1. 

Number  of  fractional  bits  for  function  values  is  8. 

Mp+ :  the  function  is  an  affine  transformation  of  an  extended  8- 
restricted  M/?-monotone  increasing  function. 

Ri  =  (EVBDD)  /  (MTBDD)  x  100.  R2  =  (EVBDD)  /  (BMD)  x  100. 
Variable  orders  of  DDs  are  produced  by  the  sifting  algorithm  [26]. 


two- variable  numeric  functions  [2].  WaveRings  in  the  table 
is 

cos  (  Vx2  +  Y 2) 

WaveRings  =  - . 

sjx2  +  Y2  +  0.25 

In  the  column  labeled  “Type  of  function’’  of  Table  3,  M p 
denotes  an  extended  8-restricted  Mp-monotone  increasing 
function,  while  M/V  denotes  an  affine  transformation  of  an 
extended  8-restricted  Mp-monotone  increasing  function. 

Two-variable  numeric  functions  whose  range  changes 
smoothly  on  a  given  domain  can  be  converted  into  ex¬ 
tended  /-restricted  Mp-monotone  increasing  functions  with 
small  p.  As  shown  in  Theorem  1,  such  functions  have 
small  EVBDDs.  In  fact,  the  two-variable  numeric  func¬ 
tions  in  Table  3  are  converted  into  8-restricted  Ml  or  MS- 
monotone  increasing  functions,  and  EVBDDs  have  many 
fewer  nodes  than  MTBDDs  and  BMDs.  Since  non-terminal 
nodes  of  EVBDD  have  weighted  1 -edges,  a  non-terminal 
node  of  EVBDD  requires  larger  memory  size  than  a  non¬ 
terminal  node  of  MTBDD  and  BMD.  However,  the  increase 
due  to  the  weighted  edges  is  negligible,  because  EVBDDs 
have  many  fewer  non-terminal  nodes  than  MTBDDs  and 
BMDs  [21], 

As  shown  in  [21],  by  converting  EVBDDs  into 
EVMDDs,  we  can  often  reduce  memory  size  and  path  length 
of  decision  diagrams.  In  the  next  section,  we  present  a  de¬ 
sign  method  that  takes  advantage  of  EVMDDs. 

4.  Two- Variable  NFGs  Based  on  EVMDDs 

4. 1  Architecture  for  NFGs 

In  DDs  based  on  the  Shannon  expansion,  function  values 
can  be  obtained  by  traversing  the  DDs  from  the  root  node 
to  a  terminal  node  [11],  [12],  In  EVBDDs  and  EVMDDs, 
function  values  can  be  obtained  as  the  sum  of  the  weights 
for  traversed  edges.  Figure  3  shows  an  NFG  based  on  an 


Input  Input  Next  variable 


(a)  Architecture  for  NFGs.  (b)  Address  computation  circuit. 

Fig.  3  Architecture  for  NFGs  based  on  EVMDDs. 

EVMDD.  It  consists  of  a  memory  to  store  an  EVMDD,  an 
address  computation  circuit  to  traverse  the  EVMDD,  and  an 
accumulator  to  compute  the  sum  of  edge  weights.  Note  that, 
for  readability,  registers,  circuits  for  initialization,  and  some 
signals  are  omitted  from  Fig.  3. 

In  Fig.  3  (a),  the  block  labeled  “Memory  for  EVMDD” 
stores  data  for  edges  in  an  EVMDD.  Data  for  an  edge  con¬ 
sists  of  a  pointer  to  the  next  node,  data  for  the  variable  of  the 
next  node,  and  an  edge  weight.  From  the  memory,  a  pointer 
to  the  next  node  and  data  for  the  next  variable  are  read  and 
fed  to  the  address  computation  circuit.  And,  an  edge  weight 
is  fed  to  the  accumulator.  The  address  computation  circuit 
produces  an  address  of  the  next  edge  from  an  address  of  the 
node  and  a  value  of  the  input  variable. 

Figure  3  (b)  shows  the  address  computation  circuit. 
Data  for  the  next  variable  consist  of  shift  data  and  mask  data. 
A  value  of  the  corresponding  input  variable  is  retrieved  by 
the  left  shifter  and  the  AND  gates.  And,  the  value  is  added 
to  the  address  of  next  node  to  generate  the  edge  address. 
This  circuit  just  selects  an  edge  to  be  traversed.  Thus,  we 
could  use  a  multiplexer.  However,  it  is  inefficient  because, 
in  EVMDDs,  each  input  variable  can  have  a  different  do¬ 
main.  Note  that  EVBDDs  require  neither  the  AND  gates 
nor  the  adder  because  in  EVBDDs,  all  input  variables  are 
binary. 

Since  the  circuit  shown  in  Fig.  3  just  traverses  an 
EVMDD  and  computes  the  sum  of  edge  weights,  it  can  eval¬ 
uate  both  one-  and  two-variable  functions  with  the  same  ar¬ 
chitecture. 

4.2  Design  Method  for  NFGs  Using  EVMDDs 

For  a  given  numeric  function,  its  domain,  and  precision,  we 
can  systematically  design  the  circuit  in  Fig.  3.  First,  con¬ 
vert  a  given  numeric  function  into  an  n-bit  precision  integer¬ 
valued  function.  Next  represent  the  integer-valued  function 
using  an  EVMDD,  and  finally  generate  HDL  code  for  the 
circuit  in  Fig.  3  from  the  EVMDD.  Since  our  NFG  directly 
realizes  the  function  table,  it  is  more  accurate  than  existing 
NFGs  using  polynomial  approximation  [7],  [  16],  [25],  [30], 
[31]. 

To  generate  memory  data  for  an  NFG  in  Fig.  3,  we  first 
assign  an  address  to  each  edge  in  an  EVMDD.  For  each  non- 
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terminal  node,  we  assign  addresses  to  edges  in  ascending  or¬ 
der  from  0-edges.  Thus,  addresses  assigned  to  0-edges  cor¬ 
respond  to  addresses  of  non-terminal  nodes.  Next,  we  gen¬ 
erate  shift  data  and  mask  data  of  each  edge.  For  an  EVMDD 
for  f(Xu,Xu- 1, . . .  ,X\),  shift  data  and  mask  data  of  an  edge 
are  computed  as  follows: 

U—  1 

shift  data  =  ^  kj,  mask  data  =  2ki  -  1 ,  and 

j=i 

bit  size  for  mask  =  max  (A:,), 

!<;<!( 

where  the  edge  points  to  a  node  representing  X;,  variable 
order  of  the  EVMDD  is  Xu,  Xu-\ , . . .  ,X\  from  the  root,  and 
kj  =  \Xj\  (j  =  1,2,...,  u). 

Example  6:  Table  4  shows  memory  data  and  initial  values 
for  the  NFG  produced  from  the  EVMDD  in  Fig.  2  (d).  This 
example  shows  how  to  compute  the  function  value  for  X  — 
(10)2  and  Y  =  (10)2  using  Table  4. 

First,  the  address  computation  circuit  produces  an  edge 
address  from  the  initial  values.  Since  the  initial  shift  data 
is  0,  the  bitwise  AND  between  the  most  significant  3  bits  of 

Table  4  Memory  data  for  the  NFG  for  2-bit  precision  2-D  norm  func¬ 
tion. 


Edge 

Shift 

Mask  data 

Address  of 

Edge 

address 

data 

(binary) 

next  node 

weight 

0 

i 

001 

8 

0 

1 

i 

001 

8 

2 

2 

0 

000 

0 

1 

3 

1 

001 

8 

2 

4 

0 

000 

0 

2 

5 

1 

001 

8 

3 

6 

0 

000 

0 

3 

7 

0 

000 

0 

4 

8 

0 

000 

0 

0 

9 

0 

000 

0 

1 

Initial  values 

Shift  data:  0  Mask  data:  111 

Address  of  node:  0  Edge  weight:  0 


input  variable  (x\  Xq  yi  )i  =  (101)2  and  the  initial  mask  data 
(111)2  is  computed  to  produce  (101)2.  Adding  the  initial 
address  of  node  0  to  the  result  of  bitwise  AND  yields  the 
first  edge  address  (101)2  =  5,  which  corresponds  to  the  edge 
5  of  the  root  node  in  Fig.  2  (d). 

Next,  data  for  address  5  is  read  from  the  memory  and 
fed  to  the  address  computation  circuit  and  the  accumula¬ 
tor.  The  accumulator  obtains  the  sum  of  the  edge  weight  3 
given  by  the  memory  and  the  initial  edge  weight  0.  In  the 
address  computation  circuit,  the  value  of  input  variable  is 
(xo  y i  yofe  =  (010)2,  which  is  shifted  1-bit  left,  and  the  bit¬ 
wise  AND  is  performed  with  the  mask  data  (001)2.  Adding 
the  result  of  bitwise  AND  0  to  the  address  of  next  node  8 
yields  the  second  edge  address  8,  which  corresponds  to  the 
edge  0  of  the  node  for  X\  in  Fig.  2  (d). 

Since  the  mask  data  0  at  the  address  8  means  arrival 
at  the  terminal  node,  adding  the  edge  weight  0  to  the  pre¬ 
vious  sum  of  edge  weights  3  yields  the  function  value  3. 

(End  of  Example) 

Memory  size  and  delay  time  of  our  NFG  depend 
largely  on  memory  size  and  path  length  of  EVMDD.  There¬ 
fore,  the  memory  minimization  algorithm  and  the  APL  min¬ 
imization  algorithm  for  MDDs  [19],  [20]  are  useful  for  pro¬ 
ducing  fast  and  compact  NFGs. 

4.3  FPGA  Implementation  Results 

We  implemented  the  NFGs  in  Fig.  3  for  the  two- 
variable  functions  shown  in  Table  3  on  an  Altera  FPGA 
(EP1S25F672C8).  To  show  the  effectiveness  of  MDDs, 
NFGs  for  EVBDDs  and  EVMDDs  were  compared.  Table  5 
shows  the  results. 

EVBDD-based  NFGs  can  achieve  higher  operating  fre¬ 
quency  than  EVMDD-based  NFGs,  because  the  EVBDD- 
based  NFGs  require  neither  AND  gates  nor  an  adder  for 
the  address  computation  circuit.  However,  as  for  the  time 
to  obtain  the  function  value  (i.e.,  latency),  EVMDD-based 


Table  5  FPGA  implementation  of  8-bit  precision  two-variable  numeric  functions. 


FPGA  device: 

Logic  synthesis  tool: 

Altera  Stratix  EP1S25F672C8 

Altera  QuartusII  7.1  (speed  optimization,  timing  requirement  of  100  MHz) 

Numeric 

functions 

EVBDD 

EVMDD 

|  Ratio  [%] 

LE 

Mem 

[bits] 

Freq. 

[MHz] 

LPL 

Delay 

[nsec] 

LE 

Mem 

[bits] 

Freq. 

[MHz] 

LPL 

Delay 

[nsec] 

Mem 

Delay 

Vx2  +  Y2 

367 

155,736 

80 

16 

201 

96 

103,080 

66 

5 

76 

66 

38 

arctan(  yVj ) 

256 

132,174 

79 

16 

204 

94 

88,760 

66 

6 

91 

67 

44 

ln(X+  l)sin(F) 

169 

142,640 

79 

16 

202 

91 

86,256 

66 

5 

76 

60 

37 

VVsin(F) 

317 

157,080 

81 

16 

198 

91 

91,404 

67 

5 

74 

58 

38 

sin(  Vv2  +  Y2) 

204 

184,968 

69 

16 

231 

93 

101,916 

65 

5 

77 

55 

33 

sin  (XY) 

365 

151,520 

81 

16 

197 

91 

90,828 

67 

5 

74 

60 

38 

X/(Y+  1) 

338 

145,782 

82 

16 

195 

91 

88,236 

67 

5 

75 

61 

38 

XY/  Vx2  +  Y2 

269 

145,200 

80 

16 

200 

94 

99,826 

67 

5 

75 

69 

37 

WaveRings 

114 

235,934 

68 

16 

234 

101 

144,892 

65 

5 

77 

61 

33 

Average 

267 

161,226 

78 

16 

207 

94 

99,466 

66 

5 

77 

62 

37 

LE:  Number  of  logic  elements  Mem:  Memory  size  Freq.:  Operating  frequency 

LPL:  Longest  path  length  of  DDs  Delay:  maximum  delay  time  =  LPL  /  Freq. 

Ratio  for  Mem  =  Mem  for  EVMDD  /  Mem  for  EVBDDxlOO 
Ratio  for  Delay  =  Delay  for  EVMDD  /  Delay  for  EVBDDx  100 
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Table  6  Performance  comparison  with  CORDIC  for  16-bit  sin(X). 


FPGA:  Xilinx  Virtex-II  XC2V 1000-6 

Logic  synthesis  tool:  Synplify  Premier  Ver.  8.5 

NFGs 

Freq.  [MHz] 

Delay  [nsec] 

CORDIC  (RTL)  [1] 

102 

157 

CORDIC  (structural)  [1] 

222 

72 

EVBDD-based 

193 

83 

EVMDD-based 

141 

43 

CORDIC  (RTL):  The  circuit  is  synthesized  from  RTL. 
CORDIC  (structural):  The  circuit  is  manually  imple¬ 
mented  to  make  it  suitable  for  the  FPGA  structure. 


NFGs  are  shorter  than  EVBDD-based  NFGs,  because  the 
path  length  of  EVMDD  is  shorter  than  that  of  EVBDD. 
The  maximum  delay  time  needed  to  obtain  a  function  value 
for  EVMDD-based  NFGs  is,  on  the  average,  only  37%  of 
that  for  EVBDD-based  NFGs.  And,  EVMDD-based  NFGs 
require,  on  the  average,  only  62%  of  the  memory  size 
needed  for  EVBDD-based  NFGs.  These  results  show  that 
EVMDDs  produce  fast  and  compact  two-variable  NFGs. 

Two-variable  functions  are  often  designed  using  a  com¬ 
bination  of  one-variable  NFGs,  multipliers,  and  adders.  For 
example,  the  compound  function  sin(  Vx2  +  Y 2)  can  be  de¬ 
signed  using  two  circuits  realizing  a 2 ,  an  adder,  a  square  root 
circuit,  and  a  sine  function  circuit.  But,  it  can  produce  a  slow 
implementation  due  to  long  path  delays.  For  such  a  com¬ 
pound  function,  two-variable  NFGs  that  can  directly  realize 
two-variable  functions  are  much  faster  [24],  To  compare  the 
EVMDD-based  NFGs  with  the  previous  two-variable  NFGs 
based  on  polynomial  approximation  [24],  we  implemented 
the  EVMDD-based  NFGs  for  Vx2  +  Y2  and  XT/  Vx2  +  Y: 
using  the  same  Altera  FPGA  (EP3SL340F1517C2)  and 
QuartusII  9.0. 

For  both  functions,  the  maximum  delay  times  of 
the  EVMDD-based  NFGs  are  21  nsec.,  while  those  of 
the  polynomial-based  two-variable  NFGs  [24]  are  35  and 
37  nsec.,  respectively.  Since  the  EVMDD-based  NFGs  are 
faster  than  the  polynomial-based  NFGs,  and  they  can  di¬ 
rectly  realize  two-variable  functions,  the  proposed  two- 
variable  NFGs  are  faster  than  the  NFGs  designed  by  any 
other  existing  approaches. 

As  mentioned  previously,  our  NFGs  can  also  realize 
one-variable  functions.  To  show  the  effectiveness  of  our 
NFGs  for  one-variable  functions,  we  compare  our  NFGs 
with  a  CORDIC  shown  in  [1],  which  is  well  known  as  a 
standard  one-variable  NFG  for  FPGA  implementation,  in 
terms  of  performance.  Table  6  shows  the  results.  From  this 
table,  we  can  see  that  the  EVMDD-based  NFG  that  is  au¬ 
tomatically  generated  from  the  function  table  evaluates  the 
sine  function  faster  than  a  manually  implemented  CORDIC. 
And,  the  EVBDD-based  NFG  is  faster  than  the  CORDIC 
synthesized  from  RTL. 

In  this  way,  we  can  implement  fast  and  compact  one- 
variable  and  two-variable  NFGs  by  using  EVMDDs. 

5.  Conclusion  and  Comments 

This  paper  introduces  a  new  class  of  integer-valued  func¬ 


tions,  called  an  /-restricted  M/j-monotone  increasing  func¬ 
tion.  It  also  derives  an  upper  bound  on  the  number  of  nodes 
in  an  EVBDD  to  represent  the  function.  EVBDDs  repre¬ 
sent  /-restricted  Mp-monotone  increasing  functions  or  their 
affine  transformations  more  compactly  than  MTBDDs  and 
BMDs  when  p  is  small. 

This  paper  also  presents  a  design  method  for  NFGs 
based  on  EVMDDs.  With  EVMDDs,  we  can  design  accu¬ 
rate,  fast,  and  compact  NFGs  for  one-  and  two-variable  nu¬ 
meric  functions.  In  FPGA  implementations  for  two-variable 
numeric  functions,  we  show  that  EVMDD-based  NFGs  re¬ 
quire,  on  the  average,  only  37%  of  the  delay  time  and  62% 
of  the  memory  size  needed  for  EVBDD-based  NFGs. 
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Appendix 

Lemma  A:  [21]  The  number  of  distinct  /2-bit  precision 
M //-monotone  increasing  functions  is  (p  +  1  )2'  1 . 

Definition  A:  A  shared  EVBDD  (SEVBDD)  is  an  exten¬ 
sion  of  an  EVBDD,  and  it  has  multiple  root  nodes  to  rep¬ 
resent  multiple  integer-valued  functions.  The  SEVBDD  is 
obtained  by  sharing  equivalent  sub-graphs  in  EVBDDs  for 
the  integer- valued  functions. 

Lemma  B:  [21]  Let  //(/,  p)  be  the  number  of  non-terminal 
nodes  in  the  SEVBDD  representing  all  the  /-bit  precision 
M //-monotone  increasing  functions,  where  the  variable  or¬ 
der  of  the  SEVBDD  is  .17  [ ,  .17... 2, . . . ,  Xo  (from  the  root  nodes 
to  the  terminal  node).  Then, 

1 

ii(Up)  =  Yj(p  +  1)2'  1 
1=1 

Proof  for  Theorem  1:  Suppose  that  an  EVBDD  for  f(X) 
is  partitioned  into  two  parts:  the  upper  and  the  lower  parts 
as  shown  in  Fig.  A-  1.  In  this  case,  the  lower  part  represents 
/-bit  precision  M //-monotone  increasing  functions,  and  the 
upper  part  represents  the  selector  function.  The  upper  part 
has  the  maximum  number  of  nodes  when  it  forms  a  com¬ 
plete  binary  tree.  That  is,  the  maximum  number  of  nodes  in 
the  upper  part  is  2"_/  -  1 . 

The  lower  part  has  the  maximum  number  of  nodes 
when  it  represents  all  the  /-bit  precision  M //-monotone  in¬ 
creasing  functions.  From  Lemma  B,  the  maximum  number 
of  nodes  in  the  lower  part  is  //(/,  p). 

Therefore,  the  number  of  non-terminal  nodes  in  the 
EVBDD  for  f(X)  is  at  most 


Fig.  A- 1  Partition  of  EVBDD. 
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l 

2"-'  -  1  +  r/(l,  p)  =  2"_/  +  Yjp  +  I  )2'”'  -  /  -  1. 

;=1 

By  adding  1  to  account  for  the  terminal  node  to  this,  we 
have  (1).  From  Lemma  A,  the  number  of  M/>monotone 
increasing  functions  that  can  be  represented  in  the  lower  part 
is  (p  +  I  )2  1 .  It  does  not  exceed  the  number  of  functions 
which  can  be  selected  by  the  upper  part:  2"  1 .  Therefore,  we 
have  the  relation:  (p  +  l)“_l  <  2"~/.  I 

Proof  for  Lemma  2:  In  the  EVBDD  for  fiX),  adding  con¬ 
stant  values  to  the  weights  of  edges  pointing  to  the  corre¬ 
sponding  sub-functions  can  produce  the  EVBDD  for  g(X). 
This  conversion  of  EVBDDs  does  not  change  the  number  of 
nodes.  I 

Proof  for  Lemma  3:  For  f(X,  Y),  any  sub-function  with 
respect  to  an  assignment  to  X  can  be  represented  by  h(Y)+b , 
where  b  is  an  integer.  Thus,  from  Definition  12,  the  lemma 
holds.  I 

Proof  for  Lemma  4: 

f(X,  Y)  =  g(X)  -  h(Y)  =  -(-g(X)  +  h(Y )) 

From  Lemma  3,  (-g(X)  +  h(Y))  is  an  extended  //-restricted 
M //-monotone  increasing  function.  Therefore,  we  have  the 

lemma.  I 

Proof  for  Lemma  5:  For  f(X,  Y),  any  sub-function  with 

respect  to  an  assignment  to  X  can  be  represented  by  a  ■  h(Y), 
where  a  is  a  real  number  satisfying  0  <  a  <  1.  Since  a  ■ 
h( 0)  =  0  and  0  <  a  •  h(Y  +  1)  -  a  ■  h(Y)  <  p  hold,  from 
Definition  11,  the  lemma  holds.  I 

Proof  for  Lemma  6:  Let  h(Y)  =  a  ■  ho(Y)  +  b.  Then, 

f(X,  Y)  =  g(X)(a  •  h0(Y)  +  b) 

=  a[g{X)-h<){Y)+h--g{X)\. 

From  Lemmas  3  and  5,  g(X)-ho(Y)+  ^  -g(X)  is  an  extended  n- 
restricted  M //-monotone  increasing  function.  Thus,  we  have 


Shinobu  Nagayama  received  the  B.S. 

and  M.E.  degrees  from  the  Meiji  University, 
Kanagawa,  Japan,  in  2000  and  2002,  respec¬ 
tively,  and  the  Ph.D.  degree  in  computer  science 
from  the  Kyushu  Institute  of  Technology,  Iizuka, 
Japan,  in  2004.  He  is  now  an  Associate  Profes¬ 
sor  at  the  Hiroshima  City  University,  Hiroshima, 
Japan.  He  received  the  Outstanding  Contribu¬ 
tion  Paper  Award  from  the  IEEE  Computer  So¬ 
ciety  Technical  Committee  on  Multiple- Valued 
Logic  (MVL-TC)  in  2005  for  a  paper  presented 
at  the  International  Symposium  on  Multiple- Valued  Logic  in  2004,  and  the 
Excellent  Paper  Award  from  the  Information  Processing  Society  of  Japan 
(IPS)  in  2006.  His  research  interest  includes  numeric  function  generators, 
decision  diagrams,  software  synthesis,  and  embedded  systems. 


Tsutomu  Sasao  received  the  B.E.,  M.E., 
and  Ph.D.  degrees  in  electronics  engineering 
from  Osaka  University,  Osaka,  Japan,  in  1972, 
1974,  and  1977,  respectively.  He  has  held 
faculty/research  positions  at  Osaka  University, 
Japan,  the  IBM  T.J.  Watson  Research  Center, 
Yorktown  Heights,  New  York,  and  the  Naval 
Postgraduate  School,  Monterey,  California.  He 
is  now  a  Professor  of  the  Department  of  Com¬ 
puter  Science  and  Electronics  at  the  Kyushu  In¬ 
stitute  of  Technology,  Iizuka,  Japan.  His  re¬ 
search  areas  include  logic  design  and  switching  theory,  representations  of 
logic  functions,  and  multiple-valued  logic.  He  has  published  more  than 
nine  books  on  logic  design,  including  Logic  Synthesis  and  Optimization, 
Representation  of  Discrete  Functions,  Switching  Theory  for  Logic  Syn¬ 
thesis,  and  Logic  Synthesis  and  Verification,  Kluwer  Academic  Publish¬ 
ers,  1993,  1996,  1999,  and  2001,  respectively.  He  has  served  as  Program 
Chairman  for  the  IEEE  International  Symposium  on  Multiple- Valued  Logic 
(ISMVL)  many  times.  Also,  he  was  the  Symposium  Chairman  of  the  28th 
ISMVL  held  in  Fukuoka,  Japan,  in  1998.  He  received  the  NIWA  Memorial 
Award  in  1979,  Distinctive  Contribution  Awards  from  the  IEEE  Computer 
Society  MVL-TC  for  papers  presented  at  ISMVLs  in  1986,  1996,  2003  and 
2004,  and  Takeda  Techno-Entrepreneurship  Award  in  2001.  He  has  served 
as  an  Associate  Editor  of  the  IEEE  Transactions  on  Computers.  He  is  a 


Jon  T.  Butler  received  the  B.E.E.  and 
M.Engr.  degrees  from  Rensselaer  Polytechnic 
Institute,  Troy,  New  York,  in  1966  and  1967,  re¬ 
spectively.  He  received  the  Ph.D.  degree  from 
The  Ohio  State  University,  Columbus,  in  1973. 
Since  1987,  he  has  been  a  professor  at  the 
Naval  Postgraduate  School,  Monterey,  Califor¬ 
nia.  From  1974  to  1987,  he  was  at  Northwestern 
University,  Evanston,  Illinois.  During  that  time, 
he  served  two  periods  of  leave  at  the  Naval  Post¬ 
graduate  School,  first  as  a  National  Research 
Council  Senior  Postdoctoral  Associate  (1980-1981)  and  second  as  the 
NAVALEX  Chair  Professor  (1985-1987).  He  served  one  period  of  leave 
as  a  foreign  visiting  professor  at  the  Kyushu  Institute  of  Technology,  Ii¬ 
zuka,  Japan.  His  research  interests  include  logic  optimization  and  multiple¬ 
valued  logic.  He  has  served  on  the  editorial  boards  of  the  IEEE  Transac¬ 
tions  on  Computers,  Computer,  and  IEEE  Computer  Society  Press.  He 
has  served  as  the  editor-in-chief  of  Computer  and  IEEE  Computer  Society 
Press.  He  received  the  Award  of  Excellence,  the  Outstanding  Contributed 
Paper  Award,  and  a  Distinctive  Contributed  Paper  Award  for  papers  pre¬ 
sented  at  the  International  Symposium  on  Multiple- Valued  Logic.  He  re¬ 
ceived  the  Distinguished  Service  Award,  two  Meritorious  Awards,  and  nine 
Certificates  of  Appreciation  for  service  to  the  IEEE  Computer  Society.  He 
is  a  fellow  of  the  IEEE. 


